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^ - Abstract 

■ The breakup of the radioactive ^^F nucleus into a proton and ^^O is studied for the reaction 

! I'^F + 208pb — > p + + 2°^Pb at 65 MeV/nucleon. The possibihty of using this reaction as 

I a test case for studying dynamical Coulomb reacceleration effects is assessed. It is shown that the 

ff^ ' reaction is dominated by elastic nuclear breakup (diffraction dissociation). 
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I. INTRODUCTION 



The Coulomb dissociation method |T[] is now a standard experimental tool as a source 
of information on radioactive capture processes of astrophysical interest. It can be shown 
that the breakup cross sections of a projectile nucleus in the Coulomb field of a target 
is proportional to the cross section for photo-dissociation [0. The radiative capture cross 
sections can be obtained via detailed balance [|T| . Thus, by measuring the Coulomb breakup 
of nuclear projectiles, specially at high bombarding energies (~ 50 — 200 MeV/nucleon), one 
can obtain information on the radiative capture cross section of interest. This goal has indeed 
been achieved in numerous experiments for the study of ^He(d, 7)^Li ^^C(n, 7)^^C 
^-^C(p, 7)^^N 0, and ^^C(a, 7)^^0 0. More recently, this method has been used exhaustively 
in the study of the reaction^Be(p, 7)^B ||, |10|, [IT], [13|, considered the most important 
one for the standard solar model fl^ . 

Under some circumstances the relation between the Coulomb breakup measurements and 
the radiative capture cross sections of interest is not so straightforward. First the radiative 
capture cross sections contain contributions of different electric and magnetic multipolarities 
which enter with different weights in the Coulomb breakup cross sections. For example, 
while the radiative E2 capture is very small in the reaction '^Be(p,7)^B within the solar 
environment, its contribution is amplified in the Coulomb breakup experiments. Separation 
of the two contributions one depends on the structure model used for the ^B nucleus |T5|, 
p!6| , p!7| , [T8[ . Secondly, the nuclear contribution to the breakup has to be separated from the 
Coulomb breakup |T^. Finally, Coulomb reacceleration effects []T^, ^ ^ have to 
be controlled. 

The Coulomb reacceleration effects are indeed one of the main concerns in extracting the 
astrophysical S-factor for the reaction '^Be(p, 7)'^B. These effects are filtered from the data 



by comparing them with dynamical calculations of Coulomb breakup M, 20, 21, 22, 23 



Due to its relevance for the application of the Coulomb dissociation method, it is desirable to 
study a system where the astrophysical S-f actors have been measured directly at the stellar 
energies and where dynamical effects can be tested. Apparently, the breakup of ^^F is a 
good candidate. The ground state of ^''F is loosely-bound (600 keV of separation energy 
into proton + ^^O) and its only excited state is one of the best halo states known sofar in 
nuclear physics, being bound by only 100 keV. Besides, the S-factor for the radiative capture 
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reaction ^^O (p, 7)^^ F has been measured down to 200 keV |2^. Indeed, the breakup of the 
weakly bound ^^F well above the barrier has been considered as an important test of the 



dynamical breakup mechanism ||2^. First theoretical analysis of this reaction has been done 
in ref. [p^ ]. 

In this work we study the breakup of ^''F into proton + ^^O, which would be relevant 
for the purpose of testing Coulomb dynamical effects by a comparison with experimental 
results. Thus, we do not consider the stripping of the proton from ^^F, or the breakup 
into other channels. We will only consider two breakup mechanisms: the Coulomb and the 
elastic nuclear breakup (or diffraction dissociation). If the proton and the ^^O fragments 
are measured in coincidence, these are the only two mechanisms of interest. We restrict our 
calculation to first order perturbation theory, as we want to learn about the feasibility of 
such experiments. We also restrict ourselves to bombarding energies of 65 MeV/nucleon, 
typical of laboratories like GANIL/France, GSI/Germany, NSCL/USA and RIKEN/ Japan, 
where previous Coulomb dissociation experiments have been carried out. 

We show that the Coulomb breakup cross sections are too small to allow for a reliable 
experimental counting rate. Also, the reaction is dominated by the elastic nuclear breakup 
mechanism which is more model dependent than the Coulomb breakup. The details of our 
calculations are presented in section 2, where we show that a single particle model is able to 
explain the astrophysical S-factor for ^^O (p,7)^^F at the lowest energies. In section 3 we 
present our calculations for Coulomb breakup and in section 4 for elastic nuclear breakup. 
Our conclusions are presented in section 5. 



II. POTENTIAL MODEL FOR ^^F 

In the present calculation, we treat ^^F as the combination of a proton and an inert ^^O 
core with spin 0. The proton-oxygen wavef unction is given by 

"^ij (r) = Yl {^lmms\jm)Yi^, (?) Xm. , (1) 

m,m,s 

where Xm^ is the spinor wavefunction. 

The radial wave functions Rij{r) are solutions of the radial Schrodinger equation for the 
proton-core motion with the potential 

-1 

+ Vc{r) (2) 



V{r) = Vo 



r ar 



1 -|- exp 



r-R 
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FIG. 1: S-factor for the reaction ^^0(p, 7)^''F. The upper (lower) data represent the capture to the 
1/2^-excited (5/2+-ground) state in ^'^F. The soHd curves are the result of a calculation for the 



direct capture with a potential model described in the text. The data are from ref. |24]. 
where 

ydf) = for r > Rc 

r 

We use a = 0.6 fm, Rc = R = 3.27 fm, Fs.o. = 0.341 fm, and tq = 1.25 fm. For the 
I"*" d-wave ground state, we use the potential depth Vq = —49.66 MeV, which reproduces 
the separation energy of 0.6 MeV. For the s-wave, we use Vq = —50.65 MeV, which 
reproduces the separation energy of 0.1 MeV of the only excited state in ^''F. 

The continuum wavefunctions are calculated with the same potential model parameters 
as the I"*" ground state. They are normalized so as to satisfy the relation 

what means, in practice, that the continuum wavefunctions UEij{r) are normalized to 
\j2mhc/T^h^k sin(A;r + 5ij) at large r, where k is the relative momentum of the fragments h 
and c (oxygen and proton, respectively). 




FIG. 2: The upper (lower) figure shows the photo-absorption cross section of ^^F for transitions 
from the ground (first excited) state to the continuum. The curves are calculated within the 
direct capture model explained in the text and are given as a function of the photon energy 
E-y = Ex = Ecm + Sp, where E^ is the excitation energy and Sp is the separation energy. 



The S-factor for the direct capture from a continuum state to the bound state, with spin 
jo, is given by 



2A[(2A + 1)!!]' k \hc 

where A = 1, or 2, for El, or E2 transitions, respectively. The electromagnetic matrix 
element Ox{E^^; Ij; /qJo) is given by 



Ox{Ecm;lj;lojo) 



ex 



47r 



/+/0+A 



X J UEij{r) ui^j^{r) r^dr, 



(6) 
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with the notation j = y/2j + 1, and ex = Zi,e{—Ac/ AaY + Zc,e{Ai,/ AaY where a = ^'^F, 
b = ^^O and c =p, respectively. 

In a similar fashion, the photo-absorption cross section for the reaction 7 + a — > b + c 
is given by 



while the photo-absorption cross section for the ground state to excited state transition 
7 a (5/2+) — > a (1/2+) is given by 



^ A IPA + 1)!!!' W |C'A(;iji;;ojo)rMB/--B.-£J, (8) 

where |CA(^iji; ^ojo) I is given by eq. with the wavefunction of the 1/2+ state replacing 
UEijir), and with Ej - Ei = 0.5 MeV. 

In figure ^ we show the S-factor for the radiative capture reaction ^^O (p, 7)^^F. The 



data points are from ref. [^. We only show the low-energy part of the spectrum, up to 
Ecm = 2 MeV. The solid curves are the result of calculations following the direct capture 
model as described above. The data are reasonably well described by the model, but the 
model overestimates the capture cross sections into the 1/2+ state for energies greater than 
1 MeV. 

In figure |^, we show the photo-absorption cross sections for transitions from the ground 
state (upper figure) and from the first excited state (lower figure) to the continuum. The 
curves are calculated within the direct capture model and are given as a function of the 
photon energy E^ = E^ = Ecm + Sp, where E^ is the excitation energy and 5*^ is the 
separation energy. Only shown are the El transition cross sections as the E2 cross sections 
are at least a factor of 10'^ smaller. The p-waves dominate the transitions for lower energies 
from the ground state to the continuum, but the f-wave contributions dominate at higher 
energies, of 2 MeV and above. The photo-absorption cross section for the transitions from 
the 1/2+ to the continuum is shown in the lower panel. 

There is also a contribution for the excitation of the 1/2+ bound-state from the ground 
state. This contribution was calculated by using eq. §. The cross section value integrated 
over the line is J a^^"^ (E^) dE^ = 7 x 10~® mb MeV. The result can be translated 

into a B(E2)-value for the down transition B{E2; 1/2+ — > 5/2+) = 61.5 fm^, which is 
rather close to the experimental value of 66.4 ± 1.4 e^ fm^ p^. This is indeed verified in 
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FIG. 3: Coulomb breakup cross section da/dE (in mb/MeV) for the reaction ^'^F (65 MeV/nucleon) 
+Pb — > ^^O + p + Pb, as a function of the proton-^^0 relative energy, in MeV. 



more elaborate calculations the reason being that the loosely-bound states in ^''F have 
a small overlap with the core, thus leading to a small core polarization. Single particle states 
are thus a good approximation for these states. 

One observes from figure ^ that the photo-absorption cross sections have maxima at 
rather high energies, (-Emax ~ 10 MeV for the g.s. and -Emax ~ 3 MeV for the excited state) 



in contrast to the photo-absorption cross sections of neutron halo nuclei [^. As explained in 



ref. p9[, the photo dissociation cross section peaks at an energy of about twice the binding 
energy of neutron halo states. The cause of the peaks moving towards the higher energies is 
the Coulomb barrier which produces a reduction of the overlap integral in eq. ^. The result 
is an effective separation energy that is much larger than in neutron halo systems with the 
same separation energies. 



III. COULOMB BREAKUP 



Since there are no data for the elastic scattering of ^''F on Pb targets at high bombarding 
energies, we construct an optical potential using an effective interaction of the M3Y type 
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FIG. 4: Coordinates used in the text to describe the F breakup. 



31| , |32| modified so as to reproduce the energy dependence of total reaction cross sections, 
i.e. m, 

t{E, S) = -I— (TNNiElab) [1 - ^a{Ei,b)] t{s) , (9) 

where to = 421 MeV fm^ is the volume integral of the M3Y interaction t{s), s is the 
nucleon-nucleon separation distance, v is the projectile velocity, a^N is the nucleon-nucleon 
cross section, and a is the real-to-imaginary ratio of the forward nucleon-nucleon scattering 
amplitude. At 65 MeV/nucleon, we use a^N = 82 fm^ and a = 0.96. 
The optical potential is given by 



U {Eiab, = j d^^i d^^i Pp (i"i)Pt (r2) t{Eiab, s) 



(10) 



where s = R + r2 — ri, and p^, (p^) is the ground state density of the target (projectile). 
For the proton, we use a gaussian density characterized by a width of 0.7 fm. For ^^O and 
208pb we use the matter densities tabulated in ref. [Q. 

Following ref. 0], the Coulomb amplitude for El transitions is given by 



J El 



2y^i'' ^^^^y^ (^^^ {j ofnol {m - mo) \jm) Oi{E^^;li;loio) Qf^il) ^M)^ 

(11) 



where /ipT is the reduced mass of the target + projectile, v is their relative velocity, 7 
(l-vVc2)-'/', and 



db b J^{qb)K^ ( ) exp [ixib)] 



(12) 



where q = 2pcm sin(^/2), pcm is the center of mass bombarding momentum and 6 is the 
scattering angle. The function J^{K^) is the cylindrical (modified) Bessel function of order 
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fi. The functions (7^ (7) are given by 



9o 



V2 

7 



9i = -9-1 = 1 ■ 



The eikonal phase, x{b), is given by 

x{h) = 21] ln{kb) — 



1 



dz Uopt{R) 



(13) 



(14) 



where rj = ZpZxe'^/hv, hk = pcm is the center of mass momentum, and R = vP+^. The 
optical potential, Uopt, in the above equation is given by eq. |10|. 
In a similar fashion, the E2 transition amplitudes are given by 



(jmjomo) 
E2 



V 30 



) {jomo2i2\jm) 02{Ecm;lj;lojo) hf,{-f) n^{q), (15) 



where fj, = m — ttiq, and 



ho = I , hi 

7 



-/i-i = - (1 + 1/7') , h2 = h. 



1 

-I—. 

7 



(16) 



The cross section for Coulomb excitation of a state with angular momentum j and excitation 
energy is obtained by averaging (and summing) over the initial (final) angular momentum 
projections: 



c 



dVtdE^ 2jo 



+ 1 



{jmjomo) I p 
El 



(jmjomo) 
E2 



(17) 



mo, m 



The El and E2 amplitudes do not interfere in the cross section after summation over 
angular momentum projections. This can easily be seen using the orthonormality condition 
of the Clebsh-Gordan coefficients: 



^ (jomo2/i|jm) (jomol/i|jm) = 0. 



(18) 



mo, m 



Thus, eq. |r7| is just a sum of cross sections for the El and E2 excitations separately: 



d^a. 



c 



d'a 



El 



+ 



d'a 



E2 



dVtdE^ dVtdE^ dVtdE^ 
where the virtual photon numbers uei and nE2 are given by 



(19) 



riEi {E^, 
nE2 {Ex, 



Ci {uj) + ^Co (^) 



7 



Zla /c\4 



7r^ 



rC2 (^) 



r 



Ci (^) + —Co (^) 



(20) 
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FIG. 5: The overlap function 0^^\ defined by eqs. 31 and 33, as a function of R, for A = 0, 1 



and 2, and for excitation energy Ex = 1.5 MeV. 
with a = 1/137 and 

Cm (^) =27i{£j jdhh Kl exp [-2 Im x (6)] ■ (21) 

Figure 3 shows the Coulomb breakup cross section da/dE (in mb/MeV) for the reaction 
^'''F (65 MeV/nucleon) + Pb — > ^^O + p + Pb, as a function of the proton-^^0 relative 
energy, in MeV. The smaller panel on the upper right corner of the figure shows the cross 
section at energies up to Ecm = 2 MeV. One sees that the El breakup mode is dominant at 
all energies. 

Although the E2 contribution to the photo-dissociation cross section is very small it 
is amplified in the Coulomb dissociation due to the large abundance of E2 photons at the 
bombarding energy considered [Q. However, the E2 contribution to da/dE is only relevant 
at Ecm ~ 3 MeV where the cross section has a peak. We also note that this cross section 
is about a factor 10^ smaller than the Coulomb breakup cross section for the reaction *B 
+ Pb — > ^Be + p + Pb at similar bombarding energies g, |, |TU|, [n|, |T|, 0. Thus, it is 
much more likely that the breakup is dominated by the nuclear interaction in this case. We 
will demonstrate this to be indeed the case in the next section. 

The angle integrated cross sections for ground state to continuum transitions are cJsi = 
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FIG. 6: Elastic nuclear breakup cross section da/dE (in mb/MeV) for the reaction ^^F (65 
MeV/nucleon) + Pb — > ^^O + p + Pb, as a function of the proton-^^0 relative energy, in 
MeV. 

0.88 mb, aE2 = 0.054 mb and (Ttotai = 0.94 mb, respectively. This should be compared to the 
cross section for the excitation of the 1/2"'' state which is 0.33 mb. These cross sections are 
small and hard to measure with reliable accuracy at present laboratory facilities. Also, since 
the strength of the Coulomb breakup is much reduced for ^^F as compared to the breakup 
of projectiles, one also expects that higher order effects, which play a role in the breakup 
of '^B, become smaller in the case of ^''F breakup. 



IV. NUCLEAR BREAKUP 



The nuclear breakup is due to the nuclear interaction 

U (R, r) = Up (rp) + Uo (ro) - Up (R) , (22) 

where Up (rp), Uq (ro) , and Up (R) are the optical potentials for the scattering of the proton, 
oxygen, and fluorine, off the target, respectively. The vectors rp, ro, and R are their 
respective coordinates relative to the target (see figure |). The fluorine-projectile potential 
is only responsible for the center of mass scattering and does not influence the ^^F breakup. 
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FIG. 7: Elastic nuclear breakup (dashed line) cross section da/dE (in mb/MeV) for the reaction 
^'^F (65 MeV/nucleon) + Pb — > ^^O + p + Pb, as a function of the proton-^^0 relative energy, in 
MeV. The solid line includes the Coulomb breakup and the nuclear-Coulomb interference. 

The DWBA matrix element for the breakup, under the assumption of spherically sym- 
metric potentials, is given by 

Tn = (R) (r) |f/p (|R - Arl) + Uo (|R - M)\ W (R)> , (23) 

where r denotes the relative coordinate between the proton and the oxygen core in ^^F, 
(R) is the ingoing and (R) is the outgoing center of mass (cm.) scattering wave- 
function, and where /5i = 16/17 and [32 = —1/17 (see fig. ||). 

For the cm. scattering at high energies one can use the eikonal approximation which 
implies 

$(-)* (R) $(+) (R) = exp [-iq ■ R + ix{b)] , (24) 
where q = k'— k is the momentum transfer to the cm. and x{b) is the eikonal phase for the 



cm. scattering, given by eq. |I4[ Since the momentum transfer is small compared to the 
total projectile momentum, we can use q = 2pcm sin(6'/2), where 6 is the scattering angle. 

It is appropriate to expand the potentials f/j (|R — /3jr|) into multipoles to exploit the 
spherical symmetry of the states in ^''F. One has then 

Ud\K- Ar I) = J2 ^) (cos e) = J2 ^y^-^f ^) ^A,(R)r,;(^) , (25) 

A X/i 
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where 



(r, R) = di (0 (^^A2r2 + i?2_2/5,ri?e) 



(26) 



and Pa (0 the Legendre polynomials. 

Employing eq. and the angular momentum algebra, we get for the case of excitation 
into a continuum state with angular momentum quantum numbers jm 



(io"io;jm;A) 



N 



zv/47r(-l)-'"° OomoA/i|jm) 



f (iolAOlil) 



Jo 



(27) 



where /i = m — mo, the function ^ is 



e^AM (r, q) = y ci'P exp [-zq ■ R + ixip)] Ty. (r, R) r,;(R) 

/■oo 

= 47r^^F,;(q) / dR R" Ja (qR) ^a (r, P) exp MR)] 
Jo 



(28) 



and^A(r,P) = E.-^? (^,^)- 

In obtaining the last result in the equation above we have used the approximation x(6) ~ 

x(P) which is here well justified, since the breakup occurs at the contact surface between 

the two nuclei, at which point z ~ 0. Another simplifying approximation is to assume that 

the momentum transfer to the cm. is mostly in the perpendicular direction, following the 

spirit of the eikonal approximation. Thus, we make the substitution Y^^{q) — > ^A^lf )' 

eq. p8| becomes 



Gx^. (r, q)=^^(-l 



A+A* 



l+lo+X 



POO 

/ dR P2 ja (qR) J=x (r, P) exp \%x{R)\ . 
Jo 



v/47r (2A + l)(A-;u)!(A + Ai)! 



(A-/x)!!(A + ^)!! 



(29) 



The (elastic) nuclear breakup cross scattering amplitude is given by (A + = even) 

e(jomo;jm.;X) / r.\ _ I^pt rp{jomo;jm;\) 



Jn 



= (_l)-™o+(A+M)/2 {j^moXfi\jm) 

X J dRR"" jxiqR)exp[txiR)] 0^^\e,^- R-lj-^jo) 



V(A-/x)!(A + /.)! 
(A-y.)!!(A + /i)!! 



(30) 
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FIG. 8: Elastic nuclear breakup (upper curve) cross section da/dfl (in b/sr) for the reaction ^^F (65 
MeV/nucleon) + Pb — > ^^O + p + Pb, as a function of the proton + ^^O center of mass scattering 
angle, in degrees. The lower curve line represents the Coulomb breakup contribution for the same 
reaction. 



where 



0)^''{E,^;R;lj;lojo) 



and the cross section is 



JoJ_ 
A2 



^0 



dr UEij{r) ui^j^ir) Tx{r,R) , 

(31) 



da 



N 



dfldE 2 jo + 



, 1 E |/] 



{jomo;jm;X) 
N 



{0) 



mo, m, A 

2 



dR R' jA {qR) exp [ix{R)] ^ (^cm; R\ I J] kjo) 



E 



dR R^ jx {qR) exp [ix{R)] Of\E,m] R] Ij] kjo) 



(32) 



with 



and 



C (Xjjo) = 2-\l E {jomoXl^\jmy 



mo, m 
(A+//=even) 



(A-^)!(A + /x)! 
[(A -/.)!! (A + /.)!! 



|2 ' 



Cf )(£;e^; R; Ij; kjo) = VC (Xjjo) Of\E,^- R; Ij; lojo) 



(33) 
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In figure ^, we plot the real part of the overlap function O^^^ as a function of R, for 
A = 0, 1 and 2, at Ecm = 1 MeV. For A = 0, only the transition from the d| ground state 
to the d| continuum states is accounted for. For A = 1, transitions to the p|, f|, and f| are 
included, while for A = 2 the transitions to the s|, d| and d| were considered. One observes 
that the overlap function O^^^ extends to large distances R between the cm. of the two 
nuclei (especially for A = 1, 2). The low distance part of this function is not relevant, due 
to the nuclear absorption at low impact parameters (included in the factor exp [ix(i?)] of 
eq. P2p. The imaginary part of 0{^^ has a similar behavior as the real part. 

In figure ^, we plot the elastic nuclear breakup cross section as a function of energy. 



obtained by integrating eq. |32| over angles. The dashed, solid and dashed-dotted curves 
represent the contributions of the A = 0, 1, and 2 multipolaritites, respectively. One observes 
that the nuclear breakup cross sections are at least two orders of magnitude larger than the 
Coulomb cross sections, as diplayed in figure ^. The elastic nuclear breakup also leads to 
fragments with larger relative energies. We have noticed this result does not vary appreciably 
with the bombarding energy. 

In figure ^ we plot the total (solid curve) breakup cross section as a function of the 
relative energy of the proton + ^^O. The cross section includes elastic nuclear and Coulomb 
breakup and Coulomb-nuclear interference. The figure is shown for the energies of interest 
for astrophysics, up to 2 MeV. The dashed curve is for the elastic nuclear breakup only. One 
clearly sees that the Coulomb breakup accounts only for a tiny fraction of the total breakup 
cross section. 

Finally, in figure || the angular distribution of the Coulomb and nuclear breakup modes 
are compared. The distribution was obtained by integration of the double differential cross 
section, eq. ^ over energy. One observes that even at the very forward angles where the 
Coulomb excitation is by comparison strongest, the nuclear breakup dominate. 



V. CONCLUSIONS 

We have analyzed the possibility of using the reaction ^'^F (65 MeV/nucleon) + Pb — > 
^^O + p + Pb for the purpose of studying the Coulomb reacceleration effects. 

We have shown that the exclusive cross section is dominated by the elastic nuclear 
breakup. This shows that this reaction is not very useful for this purpose, even if the 
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impact parameter is selected by imposing an upper cut on the center of mass scattering 
angle. In this case, the contribution of the nuclear dissociation is reduced considerably, as 
shown in ref. [^. However, the Coulomb breakup cross section is also reduced. Since the 
total Coulomb breakup cross section only amounts to fractions of milibarns, we predict a 
very small counting rate for an experiment dedicated to this study. 
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